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Abstract 

We study possible saturation effects in the total cross-sections describing interaction 
of ultrahigh energy neutrinos with nucleons. This analysis is performed within the two 
approaches, i.e. within the Golec-Biernat Wiisthoff saturation model and within the 
scheme unifying the DGLAP and BFKL dynamics incorporating non-linear screening 
effects which follow from the Balitzki-Kovchegov equation. The structure functions in 
both approaches are constrained by HERA data. It is found that screening effects affect 
extrapolation of the neutrino-nucleon total cross-sections to ultrahigh neutrino energies 
E v and reduce their magnitude by a factor equal to about 2 at E u ~ 10 12 GeV . This 
reduction becomes amplified by nuclear shadowing in the case of the neutrino-nucleus 
cross-sections and approximate estimate of this effect is performed. 



1 Introduction 

Ultrahigh energy neutrinos are one of the components of the spectrum of particles that reach 
Earth. As they interact weakly with matter their propagation through interstellar space is 
not attenuated. This is the reason why they are unique carriers of information about distant 
objects such as GRB (Gamma Ray Bursts), AGN (Active Gallactic Nuclei) etc. which are most 
probably their sources. This information can be studied by neutrino telescopes [1]. 
Attenuation of neutrinos traversing the Earth and their detection depend upon the cross- 
sections describing the interaction of neutrinos with nucleons and atomic nuclei. Ultrahigh 



energy neutrino interactions with nucleons are sensitive upon the behaviour of the nucleon 
structure functions at extremely small values of the Bjorken parameter x and relatively large 
scales Q 2 ~ [2, 3, 4, 5]. Here, as usual x = Q 2 / (2pq), where Q 2 = —q 2 with p and q denot- 
ing the four momentum of the nucleon and four momentum transfer between the leptons in the 
inelastic neutrino - nucleon interaction respectively The values of x which can be probed can 
be several orders of magnitude smaller than those which are currently accessible at HERA [6] 
and, for instance for neutrino energies E u ~ 10 12 GeV typical values of x which contribute to the 
neutrino cross-sections can be as small as x ~ 1CT 8 . Reliable estimate of the ultrahigh energy 
neutrino cross-sections does therefore require reliable extrapolation of the structure functions 
towards the region of very small values of x and large Q 2 , i.e. far beyond the region which is 
currently accessible. Existing estimates of the neutrino cross-sections with structure functions 
constrained by HERA data are based on either DGLAP [2, 4, 5] or extended BFKL [3] linear 
evolution which neglects non-linear screening effects [7, 8, 9, 10, 11, 12, 13, 14]. Those effects 
are in general expected to slow down the increase of the parton distributions and of the cross 
sections with decreasing x and to reduce their magnitude. Possible implications of screening 
effects for the estimate of the ultrahigh energy cross-sections have recently been discussed in ref. 
[15, 16, 17, 18] with somewhat conflicting conclusions. Thus in refs.[15, 16] it has been claimed 
that the screening effects should play negligible role in the estimate of the ultrahigh energy neu- 
trino cross-sections due to the dominance of relatively large scales Q 2 ~ M^. On the contrary 
results obtained in ref. [18] seem to imply that they may be significant and, moreover, when 
combined with the BFKL dynamics can even lead to enhancement of the cross-sections. 

The purpose of this paper is to present relatively detailed and realistic estimate of the impact 
of the screening effects on the cross-sections desribing interaction of ultrahigh energy neutrino 
interactions. We perform this analysis within the two frameworks which incorporate screening 
effects, i.e. the Golec-Biernat-Wiisthoff (GBW) model [19] and the unified BFKL/DGLAP 
scheme [3, 20] supplemented by the non-linear term in the corresponding evolution equations. 
This term will be obtained from the non-linear part of the Balitzki-Kovchegov (BK) equation 
[12, 13]. In both cases the parton distributions and the resulting cross-sections will be con- 
strained by the HERA data. The content of our paper is as follows. In the next Section we 
recollect basic formulas describing the deep inelastic neutrino-nucleon scattering. In Section 3 
we discuss description of deep inelastic scattering within the dipole picture and present results 
for neutrino cross sections calculated within the GBW model. Section 4 contains formulation 
of unified BFKL/DGLAP evolution equations supplemented by nonlinear screening effects and 
results for neutrino cross sections calculated within this approach. Section 5 contains summary 
and conclusions. 
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2 Basic formulas describing the deep inelastic neutrino 
scattering. 

The deep inelastic neutrino scattering is illustrated by the diagram in Fig. 1. It can proceed 

i r 




Figure 1: Deep inelastic scattering. 

through W ± or Z° exchanges that corresponds to charged current (CC) or neutral current 
(NC) interactions respectively. The charged current interactions correspond to the processes 
v + N — > l~ + X (u + N — > l + + X) with charged leptons in the final state while the neutral 
current interactions correspond to the processes u + N —>■ u + X (u + N — > u+X). The standard 
kinematical variables describing these processes are: 

s = 2ME 



Q 2 = ~q 2 

pq {U 

V= ME (1) 
where M is the nucleon mass, E denotes the neutrino energy while p and q are the four 
momenta of the nucleon and of the exchanged boson respectively. The cross-sections describing 
the deep inelastic neutrino scattering are expressed in the following way in terms of the structure 
functions F° c ' NO (x, Q 2 ), F° c ' NO (x, Q 2 ) and Ff°> NC (x, Q 2 ): 

d 2 a^f NC G 2 F ME ( M\ 



dxdy it \Q 2 + M t 



2 



1 + (1 — y) 2 j-,cc,nc/ n 2\ U 2 t? cc,nc i 



9 -Fr^(x,Q^ - »-F^{x,Q 2 )±y(l - |) xF^ c ' NC (x,Q 2 



(2) 



where Gp is the Fermi constant and Mj denotes the mass of the charged (W^ or neutral (Z°) 
gauge boson. 

In the QCD improved parton model the structure functions F2^(x,Q 2 ) are expressed in 
terms of the (scale dependent) quark and antiquark distributions [23]. Thus for the isoscalar 
target N = we have: 



Ff c (x, Q 2 ) = x[u v (x, Q 2 ) + d v (x, Q 2 )] + 2x[u(x, Q 2 ) + d(x, Q 2 ) + s(x, Q 2 ) + c(x, Q 2 )} (3) 



F° c {x, Q 2 ) = u v (x, Q 2 ) + d v (x, Q 2 ) 



(4) 



F 2 NC (x,Q 2 ) = 



(L 2 U + L 2 + Rl + R 2 ) 



x 



{x[u v (x, Q 2 ) + d v (x, Q 2 )} + 2x[u(x, Q 2 ) + d(x, Q 2 ) + s(x, Q 2 ) + c(x, Q 2 )}} (5) 



F 3 NC (x,Q 2 ) = 



(LI + L 2 -Rl- R 2 ) 



u v (x,Q 2 ) +d v (x,Q 2 )} 



(6) 



where the chiral couplings can be expressed in terms of the Weinberg angle $w 

4 

3'' 



1 - ^sin 2 6 



w 



6 r, 

L d = — 1 H — sin 6 



w 



R„ = —-sin 9 



w 



Rd = ^sin W 



(7) 



The quantities which are relevant for the quantitative description of the penetration of 
ultrahigh energy neutrinos through Earth and their detection are the neutrino cross-sections 
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integrated over available phase space at the given neutrino energy. These integrated cross- 
section are given by: 



with y = Q 2 / (xs). In equation (8) we have introduced the minimal value Q 2 min of Q 2 in order to 
stay in the deep inelastic region. In our calculations we set Q 2 min = lGeV 2 . In the "low" energy 
region s < M 2 the integrated cross-sections rises linearly with E and in this region interaction 
with valence quarks dominates. In the high energy region the contribution of valence quarks 
saturates and the energy dependence of a vD ,N (E) is driven by the small x behaviour of the 
sea quark distributions [3]. It is this part of the cross- sect ions which will be analysed in our 
paper. 

Existing numerical estimates of the ultrahigh energy cross-sections are based upon extrap- 
olation of parton distributions towards the very small x region using linear (DGLAP and/or 
BFKL) QCD evolution equations [2, 3, 4]. At small x the dominant partons are the gluons and 
the sea quark distributions are driven by the gluons through the g — > qq transitions. The linear 
QCD evolution generates indefinite increase of gluon distributions with decreasing x that im- 
plies similar increase of the sea quark distributions and of the structure functions F 2 ' N (x, Q 2 ) 
and ' (x, Q 2 ). This increase is tamed by the non-linear screening effects which lead to 
saturation [7] - [13]. Efficient way of introducing saturation can be realised using the colour 
dipole framework in which the DIS at low x is viewed as the result of the interaction of the 
colour qq dipole which the gauge bosons fluctuate to as illustrated in Fig. 2. Very succesful 



semiphenomenolgical analysis of ep DIS at low x has been performed within this framework by 
Golec-Biernat and Wusthoff [19] and in the next Section we apply this model for the estimate 
of the saturation effects in the ultrahigh energy neutrino cross-sections. 






Figure 2: Schematic representation of the dipole picture [19]. 
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3 DIS in the dipole picture and ultrahigh energy neu- 
trino interactions 



The DIS structure functions in the dipole picture can be written in the following form [19]: 

F? c L NC (x,Q 2 ) = ^J d 2 rj\z\^(r,z,Q 2 )\ 2 a d (r,x). (9) 

In this equation r denotes the transverse size of the qq dipole, z the longitudinal momen- 
tum fraction carried by a quark and Vvif (r, z, Q 2 ) are proportional to the wave functions of 
the (virtual) charged or neutral gauge bosons corresponding to their transverse or longitudi- 
nal polarisation (F^ c,NC (x, Q 2 ) = F° c ' NC (x,Q 2 ) - F° c ' NC (x, Q 2 )). Explicit expressions for 
$tl ( r ? z i Q 2 ) are gi ven below. The cross-section a d (r, x) describes interaction of the colour qq 
dipole with the nucleon. In the GBW model a d (r,x) has the following form: 



a d (r,x) = (7 



1 — exp 



y»2 



(10) 



The most crucial element in this model is the adoption of the x-dependent saturation radius 
Ro(x) which scales the qq separation in the dipole cross-section. Saturation radius is the 
decreasing function with decreasing x and is parametrised as below: 



Qo \ x 

with Ql = lGeV 2 . The three parameters of the model cr , A and x were fitted to inclusive DIS 
data from HERA for x < 0.01. We shall use the following values a = 29.12m&, A = 0.2777 
and x = 0.41 x 10~ 4 which were obtained from the fit with four flavours. 



In the limit r — > oo we have a d — > o"o , i-e. the dipole cross-section exhibits saturation prop- 
erty. The fact that the dipole cross-section is limited by the energy independent cross-section 
can be regarded as the unitarity bound. In the limit r — > the dipole cross-section vanishes 
reflecting the colour transparency. 



The GBW model, which has proved to be phenomenologically very succesful in describing 
HERA data and has embodied saturation can be used for the estimate of the UHE neutrino 
cross-sections. In our calculation of UHE neutrino-nucleon cross-section we consider only four 
flavours (u, d, s, c). The corresponding qq dipoles which contribute to Cabibbo favoured transi- 
tions are ud(du), cs(sc) for charged currents and uu, dd, cc, ss for neutral currents respectively. 
In our calculations we shall assume massless quarks. This approximation is reasonable for very 
high energy neutrinos. Possible contrtibution of heavy quarks (b, t) where the mass parameters 
cannot be neglected is found to be relatively small [3]. 
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The dipole model describes well the deep inelastic scattering at small x, but it becomes 
inaccurate at large and moderately small values of x. This is closely linked with the fact that 
it neglects theoretical expectations concerning behaviour of quark distributions in the x — + 1 
limit. This effect can be approximately taken into account by multiplying the structure func- 
tions by a factor (1 — a;) 2 ™ 3-1 which follows from the constitutent counting rule where n s denotes 
the number of spectator quarks. Since the dipole model represents the sea quark contribution 
we set n s = 4. 

The functions (r, z, Q 2 ) corresponding to the sum over dipoles corresponding to massles 
quarks are given by the following formulas: 

\^(r,z,Q 2 )\ 2 = ^[z 2 + (1 - zf]Q 2 Kl{Qr) (12) 

\^(r,z,Q 2 )\ 2 = - 2 z\\ - z) 2 Q 2 K 2 (Qr) (13) 

\4>l(r, z, Q 2 ) | 2 = ^(L 2 + L\ + Rl + R 2 )[z 2 + (1 - zf]Q 2 Kl{Qr) (14) 

0£(r, z, Q 2 )\ 2 = - 2 (Ll + L 2 + Rl + R 2 )z 2 (l - z) 2 Q 2 K 2 (Qr) (15) 

where 

g 2 = ^(i-^)g 2 (16) 

and K 0t i(u) are the Mc Donald's functions. 

In Figures 3 and 4 we show results for the ultrahigh neutrino cross-sections calculated 
within the GBW saturation model and confront them with the estimate based upon the unified 
BFKL/DGLAP framework which ignored saturation effects [3]. We can see that the cross- 
sections calculated within the GBW model are at ultrahigh energies E ~ lCPGeV about a 
factor two smaller than those which were estimated within the scheme incorporating BFKL 
and DGLAP evolution without screening corrections. 
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Figure 3: The prediction for the neutrino nucleon CC cross section obtained from the 
GBW saturation model. For comparison we also show results based on the (linear) unified 
BFKL/DGLAP evolution. 
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It should however be remembered that the GBW model has been tested phenomenologically 
for relatively low values of Q 2 and that it requires corrections incorporating DGLAP evolution 
[21]. It is therefore necessary to perform an estimate of the neutrino cross-sections within the 
scheme that would include saturation effects together with complete QCD evolution as it will 
be described in the next Section. 

4 Unified BFKL/DGLAP evolution with nonlinear screen- 
ing effects and ultrahigh neutrino cross-sections 

We have shown in the previous Section that the ultrahigh energy neutrino cross-sections based 
upon the GBW saturation model are at very high neutrino energies (E > 10 12 GeV ) about 
a factor two smaller than those calculated from the linear QCD evolution equations. It can 
however be expected that part of this reduction may just be caused by the fact that the 
GBW model does not correctly include the QCD evolution effects [21]. Let us recall that 
the dominant contribution to the neutrino cross sections comes from the region Q 2 ~ 
where the simple GBW model may not be sufficiently accurate. It would therefore be desirable 
to discuss the ultrahigh energy cross-sections within the framework which contains both the 
QCD evolution effects and saturation. It is also of course very important that this framework 
should be based upon realistic parton distributions constrained by HERA data. The framework 
should also contain all the QCD expectations concerning the small x behaviour which follow 
from the BFKL dynamics as discussed, for instance, in ref. [20]. We shall therefore use the 
scheme developed in ref. [20] containing the unified BFKL/DGLAP dynamics with subleading 
BFKL effects taken into account and supplement it by the screening contributions for the gluon 
distributions. To be precise we shall include the non-linear screening term in the corresponding 
equation for the unintegrated gluon distribution f(x, k 2 ), where k 2 is the transverse momentum 
squared of the gluon [20]. The structure functions are then calculated from the unintegrated 
gluon distributions using the k t factorisation prescription [25]. The extended system of the 
evolution equations with screening effects included then reads: 

f(x,k 2 ) = f(°\x,k 2 ) 
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- - k2 i) 2 Iff [f (B '<*• H 2 « i7) 

The first three lines in equation (17) describe the linear unified BFKL/DGLAP evolution [20]. 
Thus the second line of this equation corresponds to the BFKL evolution [7, 26]. The constraint 
— k' 2 j reflects the so called consistency constraint [27] which generates dominant part 
of the subleading BFKL corrections [28, 29]. The two terms in the third line in eq. (17) 
correspond to the DGLAP effects generated by that part of the splitting function P gg (z) which 
is not singular in the limit z — > and by the quarks respectively, with E(x, k 2 ) corresponding 
to the singlet quark distributions 

E(x, k 2 ) = (? + <?) + c + c = V(x, k 2 ) + + S c (18) 

where V, S u d s and denote the valence, the light sea quark and the charmed quark distributions 
respectively. The inhomogeneous term f(°\x,k 2 ) is defined in terms of the input (integrated) 
gluon distribution: 

/<°>(s, k 2 ) = ^1 £ dzP 99 {z) X -g g, A; 2 ) (19) 

The nonlinear screening contribution is given by the last term in equation (17) where R 
denotes the radius within which the gluons are expected to be concentrated. The structure of 
this contribution follows from the Balitzki-Kovchegov equation [12, 13] adapted to the uninte- 
grated gluon distribution f(x, k 2 ) [24] and the details concerning the structure of this term are 
briefly discussed in the Appendix A. The sea quark distributions which describe the structure 
functions are calculated from the unintegrated gluon distributions using k t factorisation. 



Unlike the case of the leading ln(l/x) approximation equation (17) cannot be reduced to the 
evolution equation in ln(l/x) that makes its numerical analysis rather cumbersome. This com- 
plication comes from the fact that the subleading BFKL effects and the non-singular DGLAP 
contribution introduce nontrivial z dependence of the kernel. We have however observed that at 
small x linear version of equation (17) can be very well approximated by an (effective) evolution 
equation in ln(l/x) with the boundary condition provided at some moderately small value of 
x (i.e. x = xo ~ 0.01). The latter is obtained from the solution of linear version of equation 
(17) in the region x > x . To be precise we have observed that at small x we can make the 
following approximations: 

I. The consistency constraint <d(k 2 /z—k' 2 ) of the BFKL kernel responsible for the subleading 
BFKL effects is replaced by the following effective (z independent) term 

&(k 2 /z - k' 2 ) -> &(k 2 - k' 2 ) + (^-] e " &(k' 2 - k 2 ) (20) 
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This replacement is motivated by the structure of the consistency constraint in the mo- 
ment space, i.e. 



CO 



£ ^ Q ^/ z _ k' 2 ) = &(k 2 - k' 2 ) + (JLy e(k' 2 - k 2 ) 



(21) 



2. We make the following replacement: 

[ ^[zP m (z) - 2N c ]f [- z , k' 2 ) - P gg {u = 0)f(x, k' 2 ) 

where P gg {ui) is a moment function of zP gg (z) — 2N C , i.e. 

f 1 dz 

Pp fl H= / -z"[(zP gg (z)-2N c ] 
Jo z 



(22) 



(23) 



This approximation corresponds to keeping the leading term in the expansion of P gg (u) 
around uj = that is a standard approximation at low x [30]. 

3. We neglect the quark contribution in the right hand side of equation (17). 

Using these approximations equation (17) can be rearranged into the following form: 

df(x,k 2 ) 



dln(l/x) 



= K'" <8 / 



(24) 



where 



Kl" ® / = 



2tt 



k- 2 



"2 



2 k' 2 





fix, k' 2 ) 


0(A;"2 _ fc/2) + ^y e,f Q(k' 2 - k" 2 ) 


- /(*, k" 2 ) 


\k' 2 -k" 2 \ 



+ 



[4k 14 + k" 4 ]^ 



(25) 
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K s (k\ k" 2 ) = 5(k 2 - k" 2 ) + Q(k 2 - k" 2 )P gg (0)^lexp [P 99 (0)m 2 ) ~ (26) 



and 



Following ref. [20] we set k 2 = lGeV 2 . The method of the solution of equation (24) is described 
in the Appendix B. 

At first we solved the linear version of equation (24) with the non-linear term neglected 
starting from the initial conditions at x — obtained from the solution of the exact 

equation. The parameter oj e ff was then obtained by fitting the solution of the linear ver- 
sion of approximate equation (24) to the solution of linear version of exact equation (17) of 
the unified BFKL/DGLAP framework. This procedure gives u e ff = 0.2. It turns out that 
the solution of the linear version of the approximate equation (24) reproduces the solution of 
the linear version of equation (17) within 3% accuracy in the region 10~ 2 > x > 10~ 8 and 
2GeV 2 < k 2 < 10 6 GeV 2 ). 

We next solved the non-linear equation (24) setting R = AGeV -1 . The quark distributions 
defining the structure functions F^' NC were calculated from the k t factorisation [3, 20] 

irj ^S b ™(z,k 2 ,Q 2 )f(* z ,k 2 ) (28) 

where the impact factors corresponding to the quark box contributions to gluon-boson fusion 
process are the same as those used in ref. [20] (see also [31]), i.e.: 

Q 2 r 1 .**^ 



S b q ox (z, k , Q 2 ) = J df3d 2 n'a s b{z - z ) 



i \ 2 

K K — K 



r + (! - PY\ [p- - -T7- ) + H - l( r n 1 - ■ i [ — - — ) ) (29) 
where k' — k — (1 — f3)k and 




D lq = k 2 + - (3)Q 2 + m\ 
D lq = (K-k) 2 + f3(l-f3)Q 2 + ml 
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-1 

(30) 

To be precise in the calculation of the (effective) quark distributions appearing in the charged 
current structure function we use the impact factors (29) corresponding to the massles quarks 
and the (charmed) quark mass effects are included in the threshold factors: 

a c ,(A: 2 )=(l + ^^)" (31) 

The k t factorisation formulas (28,29) contain subleading ln(l/x) effects coming from the exact 
kinematics of the gluon-boson fusion process [32]. Complete NLO corrections to the impact 
factors are discussed in [33]. In the impact factors corresponding to the neutral currents we use 
equation (29) with m u = = m s = and m c = lAGeV. We also include non-perturbative 
contributions according to the prescription defined in ref. [20]. The valence quark distribu- 
tions were taken from ref. [22]. In Figures 5 and 6 we show results of our calculation for 
<j cc = a v C (E) and a N c = ^^{E) with and without screening corrections included and 
confront them with our previous estimate based upon the GBW model. We can see that 



zo = 



1 + 



k' z + nig 
(3(1 -P)Q 2 



k 2 

~Q 2 
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Figure 5: The prediction for the neutrino nucleon CC cross section obtained from unified 
BFKL/DGLAP equation supplemented by screening effects. For comparison we also present 
results based on the GBW saturation model and the linear unified BFKL/DGLAP evolutions 

at ultrahigh energies the cross-sections calculated within the unified BFKL/DGLAP frame- 
work supplemented by screening effects are bigger than those calculated from the simple GBW 
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Figure 6: As for Fig. 5, but for the NC interactions. 



model. The resulting cross-sections are still appreciably smaller than the cross-sections calcu- 
lated within the linear BFKL/DGLAP framework with the screening effects neglected. 

Reduction of the magnitude of the neutrino cross-section is the consequence of the fact that 
the non-linear screening effects slow down increase of the structure functions with decreasing x. 
In Fig. 7 we show the charged current structure function F% c {x, Q 2 ) plotted as the function of 
x at Q 2 = with and without screening corrections included. We can see that the screening 
effects reduce the magnitude of F 2 c * c '(a;, Q 2 = M^) at x — 1CT 8 by almost a factor equal to two. 
The screening effects in structure functions are generated through k t factorisation by the screen- 
ing effects in the unintegrated gluon distribution f(x, k 2 ) which satisfies the non-linear equation 
(24). The non-linear screening corrections generate the critical line Q 2 {x) which increases with 
decreasing x which divides the k 2 , x plane into the two regions. In the region k 2 < Ql(x) the un- 
integrated gluon distribution saturates, i.e. f(k 2 ,x) ~ R 2 k 2 h(x,k 2 ), where h(x,k 2 ) is a slowly 
varying function of x and k 2 . In this region the unintegrated gluon distribution becomes then 
much smaller than the solution fi(k 2 ,x) of the linear version of equation (24) which behaves 
approximately as fi(k 2 ,x) ~ x~ x with A ~ 0.3 over the entire region of k 2 [3, 20]. In the region 
k 2 > Q 2 c {x) the non-linear screening contribution in the right hand side of equation (24) becomes 
less important than the linear term and can be neglected for k 2 » Q 2 (x). The magnitude of 
the unintegrated distribution continues to be significantly smaller than fi(x, k 2 ) over substan- 
tial range of k 2 . To be precise we have f(k 2 ,x) ~ f{Q 2 {x),x) with f(Q 2 (x),x) < fi(k 2 ,x) in 
the region 

Q 2 (x) >k 2 > Ql(x) (32) 
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Figure 7: The F 2 CC (x, Q 2 ) structure function obtained from the unified BFKL/DGLAP equa- 
tion supplemented by screening efects compared to results based on the linear BFKL/DGLAP 
evolution. The function F 2 cc (x, Q 2 ) is plotted as the function of x for Q 2 = Mjy. 

where Q 2 (x) = Q 4 c (x)/A 2 QCD [34]. 1 It should be noted that Q 2 (x) » Q 2 c (x). The screening 
effects do therefore significantly reduce the corresponding contribution to the k t factorisation 
integrals (28) coming from the region k 2 < k 2 < Q 2 (x). The integral over this region gives 

CC N C o 

of course part of the leading twist contribution to the structure functions F 2 ' (x, Q ) that 
does not vanish at large Q 2 . This result that the screening effects at the structure function F 2 
are appreciable even at such a large value of Q 2 ~ comes therefore from the fact that the 

CC N C o 

screening effects contribute to the leading twist part of F 2 ' (x, Q ). 

The fact that the screening effects at F 2 ,N (x, Q 2 ) can be important at Q 2 ~ and 
very small x (x ~ 1CT 8 ) implies that they may in turn have non-negligible influence on the 
ultrahigh energy neutrino cross-sections. It is this fact which makes our results significantly 
different from those presented in refs. [16, 17] where the saturation effects were confined to 
the modification of the structure functions in the saturation region Q 2 < Q 2 {x) only. The 
corresponding contribution to the UHE neutrino cross-section coming from the integral over 
this region in equation (8) is very small and so modifications of the structure functions in the 
saturation region alone have negligible impact on the UHE cross-sections [16, 17] . 

The fact that the cross-sections are sensitive upon the behaviour at very small x and large 
scales Q 2 ~ implies that the effects which are formally subleading in ln(l/x) but can signif- 
icantly affect both the ln(l/x) and Q 2 evolutions cannot be neglected. We illustrate this point 

1 Condition (32) has a simple origin. It comes from the fact that possible scaling violations in the region 
k 2 > Q 2 (x) which modify boundary condition provided along the critical line Q\ (x) are approximately controlled 
by the 'evolution length' £(fc 2 ,x) ~ a s (Q 2 (x))ln(k 2 / 'Ql(x)). Condition (32) is equivalent to the requirement 
C(fc 2 ,x) << 1 
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in Fig. 8 where we show F 2 cc, (x, Q 2 ) for Q 2 = calculated within three approximations: 
the leading ln(l/x) BFKL framework, the extended BFKL framework which includes the sub- 
leading ln(l/x) effects generated by the consistency constraint and the unified BFKL/DGLAP 
scheme which includes besides the BFKL dynamics with subleading effects also the complete 
DGLAP evolution. The latter two frameworks contain effects which are subleading at low x. 
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1 - 

le-08 le-06 0.0001 

x 

Figure 8: The comparison of the F 2 CC * structure function calculated for Q 2 = obtained 
from unified BFKL/DGLAP evolution, extended BFKL evolution and LO BFKL evolution. 
Extended BFKL evolution corresponds to the BFKL equation with subleading effects generated 
by the consistency constraint but without the DGLAP effects. 

For simplicity of presentation the non-linear screening effects are neglected in all three cases. 
We can see that both subleading ln(l/x) effects play very important role and significantly re- 
duce the magnitude of the structure function at large scale Q 2 ~ and very small values of 
x ~ 10~ 8 . Including these effects together with screening contribution is therefore important 
for getting a reliable extrapolation of the structure functions into the region of very small values 
of x and large scales. 

Discussion of the cross-sections performed so far concerned screening effects on a nucleon 
target. In the case of the neutrino- nucleus inelastic scattering further reduction of the magni- 
tude of the total neutrino cross-sections due to nuclear shadowing is expected [36]. In order to 
perform an indicative estimate of the possible nuclear shadowing effects for different values of 
the atomic numbers A we just modify the strength of the non-linear of the non-linear term in 
eq. (17) by a factor A 1 / 3 . In Fig. 9 we show our results for the normalised neutrino- nucleus 
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cross sections for different values of the atomic number A varying from A = 12 to A = 207. For 
comparison we show results for the neutrino-nucleon cross-section with and without screening 
effects. We see from this Figure that the nuclear shadowing can lead to further reduction of 
the cross-section. 
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Figure 9: The prediction for the neutrino nucleus CC cross section obtained from unified 
BFKL/DGLAP equation suplemented by screening effects. Cross section is calculated for 
diffrent atomic numbers and normalised to nucleon. For comparison we also present results for 
neutrino nucleon CC cross section based on the (linear) unified BFKL/DGLAP evolution. 



5 Summary and conclusions. 

In this paper we have performed analysis of possible implications of the screening effects on 
the extrapolation of the neutrino-nucleon cross sections towards the ultrahigh energy region. 
Behaviour of the cross-sections in this region probes the structure functions at very small val- 
ues of x and relatively large scales Q 2 ~ M^ z . The values of x which can be probed can 
be as small as 1CT 8 and it may be expected that parton densities in this ultra small x region 
should be affected by non-linear screening effects which tame the indefinite increase of parton 
distributions generated by linear (BFKL and/or DGLAP) QCD evolution. At first we have 
performed an estimate of the total neutrino-nucleon cross sections within the Golec Biernat 
- Wusthoff saturatuion model. In this model the deep inelastic lepton scattering is viewed 
as the result of the interaction of the colour qq dipoles which the gauge boson fluctuates to. 
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Important ingredient of the model is the fact that it incorporates saturation property of the 
total dipole-nucleon proton cross-section at large transverse separations between constituents 
of the dipole. We have found that the neutrino total cross-sections obtained within this model 
are at ultrahigh neutrino energies significantly smaller than those estimated from the linear 
BFKL/DGLAP evolution which neglects screening effects. We have observed however that 
part of this reduction might have been caused by the fact that the GBW model did not include 
the QCD evolution effects and so it was not sufficiently accurate at very small values of x and 
large value of the scale Q 2 ~ M^ z . In order to overcome this potential deficiency of the GBW 
model we have performed an estimate of the cross section within the more elaborate framework 
based on the unified BFKL/DGLAP scheme supplemented by the non-linear screening effects. 
Contrary to the simple GBW model this framework contained complete BFKL and DGLAP 
evolution including the subleading BFKL contributions. We have shown that all these effects 
are important in the region of very small values of x and large scales which is relevant for the 
interactions of ultrahigh energy neutrinos. The non-linear screening effects were still found 
to reduce appreciably the neutrino cross-sections at ultrahigh energies yet their effect turn 
out to be milder than in the case of the simple GBW model. We have also presented a very 
approximate estimate of the nuclear shadowing effects on further reduction of the cross-sections. 

To summarise we have shown that the screening effects may play non-negligible role in the 
extrapolation of the neutrino cross-sections towards the ultrahigh energies. 
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In this Appendix we derive the non-linear shadowing term in the right hand side of eq. (24) 
starting from the Balitzki-Kovchegov equation. The basic quantity within this framework is 
the number of the colour dipoles N(r, b, x) in a nucleon where r denotes the transverse size of 
the dipole and b is the impact parameter. The quantity N(r, b, x) is closely related to the total 
cross-section a(r,x) describing interaction of the qq colour dipole of transverse size r with a 
nucleon 



The dipole cross-section is related to the unintegrated gluon distribution f(x, k 2 ) [35]: 



Appendix A 




(33) 




(34) 



17 



For simplicity we regard a s as fixed parameter and will put its argument at the end. It is 
convenient to introduce the functions N(l, b, x) and n(l,x) defined as below: 



N(l,b,x) = 



d 2 Y 

2nr 2 



exp[i\r]N(r, b, x) 



«(u) = /^ab,») 

From equations (33,34,35, 36) we get: 

n 2 roc ^2 fc2 

where we have used the following relation: 

r°° A r u 
/ -Mlr)[l - Jo(kr)} = &(k 2 - l 2 )ln(-) 
Jo r <• 

Equation (37) implies the following local relation between h(l,x) and f(x,l 2 ): 

{l-l^)H 2 n{l,x) = ^ff{x,l 2 ) 



(35) 
(36) 

(37) 



(38) 



(39) 



In the large iV c limit the function N(r, b, x) satisfies the Balitzki-Kovchegov equation [12, 13]: 

iV(r i, b, x) = N (r 01 , b, x) 



+ 



a.N r f 1 dz 



2n 



[ -{-2ln r -^N(v 01 ,b,z) 

Jx Z p 2 



°°d 2 v 2 r 2 m 



2 A' ( r 02 , b + ir 12 , zj-N (r 02 , b + ^r 12 , zj N (r 12 , b - ^r 20 , z 



} (40) 



The term linear in A^ in the right hand side of equation (40) corresponds to the right hand 
side of the BFKL equation (in transverse coordinate space) in the leading ln{l/x) approximation 
and the nonlinear term describes screening effects. Taking the Fourier-Bessel transform of both 
sides of eq.(40), integrating over d 2 b using an approximation b » l/2r 20 and b » l/2r 10 in 
the nonlinear term and assuming the following factorisation: 



A> 2 (l,b,^) =n(l,x)S(b) 



with 



/ 



d 2 bS(b) = 1 



(41) 



(42) 
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we get 



where 



l-nil.r) l-n' ] (l-..r) ■ ' 



7T 



dz 



K®l 2 n 



7lR 2 



l 2 n 2 (l,z) 



irR 2 



d 2 bS 2 (b) 



(43) 



(44) 



The kernel K in equation (43) is the LO BFKL kernel. Using equations (37) and (39) we 
transform equation (43) into an equation for the unintegrated gluon distribution: 



where following ref. [21] we set argument of a s equal to A;' 2 in the non-linear term. Finally 
supplementing the linear evolution by the subleading BFKL effects generated by the consistency 
constraint and the DGLAP contributions resulting from the non-singular parts of the P gg 
splitting function and quarks we get equation (17). Following ref. [20] we set argument of a s 
in the linear term equal to k 2 . 



Appendix B 

In order to solve equation (24) we use the Tchebyshev interpolation formula: 

2 N ~ 1 

fc 2 ) = ^E v n T n (r t )T n (r k ,)f(x, k 2 ) (46) 

n,i=0 

where T n (z) are the Tchebyshev polynomials, vo = 1/2 and v n = 1 for n > 0. The variables 
Tfc2, Tj and k 2 are defined as 

In 

n = cos (^T^ (48) 
k 2 = k max k ( ^p-Y (49) 



ko 

where we set k^ = lGeV 2 , k 2 max = 10 6 GeV 2 . Combining the Tchebyshev interpolation formula 
(46) with equation (24) we reduce this equation to a system of the non-linear differential 
equations for the functions f(x,k 2 ). This system is solved using the standard Runge-Kutta 
method and the function f(x,k 2 ) calculated from (46) for arbitrary value k 2 in the region 
fco < k 2 < k^x- In the region k 2 > k^^ which gives negligible contribution anyway we 
approximate the function f(x, k 2 ) by f(x, k^^). 
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